Precision Determination of Heavy Quark Masses and the Strong Coupling Constant 
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We present a new QCD sum rule with high sensitivity to the continuum regions of charm and 
bottom quark pair production. Combining this sum rule with existing ones yields very stable results 
for the MS quark masses, rhc{mc) and rhi,{rhi,). We introduce a phenomenological parametrization 
of the continuum interpolating smoothly between the pseudoscalar threshold and asymptotic quark 
regions. Comparison of our approach with recent BES data allows for a robust theoretical error esti- 
mate. The parametric uncertainty due to is reduced by performing a simultaneous fit to the most 
precise sum rules and other high precision observables. This includes a new evaluation of the lifetime 
of the r lepton, Tt, serving as a strong constraint on as. Our results are rhc{Thc) = 1.2891q gls GeV, 



= 4:.207tQ°li GeV (with a correlation of 29%) 



and as{Mz)[Tr] = 0.1221 



-0.0026 
-0.0023- 



m 
> 



o . 

(N ■ 

O ' 

Oh; 

I 

Oh: 



X 



The determination of the fundamental Standard Model 
(SM) parameters is important not only in its own right, 
but also as a SM test when results from various sources 
are compared. Such comparisons can foster our under- 
standing of SM dynamics (such as strong QCD effects) 
or may ultimately lead to hints of new physics beyond 
the SM. Moreover, precise values of the SM parameters 
can be compared against the predictions of more funda- 
mental theories. Well known examples are Grand Unified 
Theories which typically predict values for the strong 
coupling constant, as, and the mass ratio mb/rriT. 

It is generally difficult to obtain reliable informa- 
tion on quark masses. The Particle Data Group 2j 
lists only ranges for their values, indicating a lack of 
confidence in the theoretical methods used to evaluate 
them. Indeed, is quite large at the mass scales of 
the bottom and charm quarks, questioning the conver- 
gence of perturbative QCD (PQCD). Furthermore, non- 
perturbative (power suppressed) effects governed by the 
scale Aqcd ~ 0.5 GeV could be large, thus compromis- 
ing reliable calculations. Two types of conditions arc 
known to improve the situation: high energy or inclu- 
siveness. As an example for the former, as and nib 
can be determined at LEP energies using PQCD. This 
yields a^iMz) = 0.1200 ± 0.0028 |3| with very little the- 
oretical uncertainty, but b quark effects are small and 
rhb{Mz) = 2.67 ± 0.50 GeV 4] is not weU constrained. 

In this Letter, we compute as from r-r, by definition an 
inclusive quantity and known to be quite insensitive to 
effects from non-perturbative QCD (NPQCD) 0. Like- 
wise, we use a set of inclusive QCD sum rules to derive 
values for rhc{rhc) and mb(fhb). One of these sum rules is 
new, and its use together with existing ones 0, Q proves 
to be a powerful tool to constrain the continuum region of 
quark pair production. This will be particularly helpful 
for the case of the b quark for which precise measurements 



of R{s) (the inclusive hadronic cross section normalized 
to the leptonic point cross section) or of Rb{s) (exclusive 
cross section for hh pairs) are unavailable. 

On the basis of an unsubtracted dispersion relation 
(UDR) it was shown in Ref. that knowledge of TOc, 
TOfc, and as is sufficient to compute the charm and bottom 
quark contributions to the QED coupling a{-Jt = Mz), a 
vital parameter entering the analysis of the very high pre- 
cision LEP 1 and SLC data. Or conversely, comparison 
of this UDR with the more traditional approach using a 
subtracted dispersion relation (SDR) offers information 
on nic and mf,. The resulting equation relates an inclu- 
sive integrated cross section to a difference of vacuum 
polarization tensors, viz. 
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Eq. JQ) defines a continuous set of sum rules parametrized 
by t, where the limit < ^ coincides with the first mo- 
ment of IIq{t). Similarly, there is a sum rule. 
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for each higher moment, Mn, as well [ESSEiOl- We 
now take the opposite limit in Eq. J^l, i ^ oo, and regu- 
larize the divergent expression (which will render A^o < 

0, 
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Qq and nq are the quark charge and the number of active 
flavors. Using expressions derived in Refs. we can 
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now explicitly write down the sum rule Q for t — s- cx): 
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Here, Af^ and are the mass and the electronic partial 
width of resonance R, and i?™"* denote the continuum 
regions integrated from M = Mq± for b and M = A/^jn 
for c. The C(2) terms arose from the regularization (PJ 
which together with the scale choices rhq = rhq{rhq) and 
ds = ds(TOq) eliminates (resums) all logarithmic terms in 
Eq. Unlike in any of the sum rules l(2Jl, i?™"' appears 
unsuppressed in Eq. Q so that rhq varies exponentially 
with the experimental information on the resonances. An 
optimal approach to compute a{Mz) would first identify 
the sum rule most sensitive to m^, and then use the value 
so obtained in theoretical expressions such as the one 
presented in Ref. . We will use Eq. to constrain the 
continuum region and work with the following ansatz: 
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where now ds = ds(2M), s' = s + 4{m'^^{2M)- AP), and 

ds(2M) 
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We will use the form in the second line (applying it to 
all moments) of Eq. jSJl with the corresponding change in 
the regularization in Eq. Q . This keeps only the leading 
logarithms resummed but allows for an analytical inte- 
gration. Eq. (jSJ coincides asymptotically with the pre- 
dictions of PQCD for massless quarks and interpolates 
smoothly between the vanishing phase space at the pseu- 
doscalar threshold and the strong onset of fermion pair 
production. The quark parton model predicts A? = 1, 
while from third order massive QCD corrections [lj| one 
expects A3 > 1 (in agreement with our results). But 
unlike when PQCD is apphed to R{s) directly and rela- 
tively close to the resonance region, we minimize the ex- 
posure to local quark-hadron duality violations by using 
QCD inclusively and by merely requiring stable results 
across the moments. No claim is being made about the 
local shape of Rq — we only need theoretical information 
about global averages. 



We use the narrow resonance data Q listed in Table U 
as the only experimental input. The wider resonances in 
the continuum region are assumed to be accounted for by 
our ansatz (jSl because (i) they decay almost exclusively 
into flavored hadrons; (ii) they interfere with the non- 
resonating part of the continuum rendering a common 
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Mr [GeV] 


rie [keV] 


R 


Mr [GeV] 


[keV] 


J/* 
*(25') 


3.09687 
3.68596 


5.26 (37) 
2.19 (15) 


T(IS') 
T(2S) 
T{3S) 


9.46030 
10.02326 
10.35520 


1.320 (50) 
0.520 (32) 
0.476 (78) 



TABLE I: Resonance data used in the analysis. The un- 
certainties from the resonance masses are negligible. 

treatment virtually impossible; (iii) the J-function ap- 
proximation (which is perfect for the narrow resonances) 
becomes successively worse; (iv) the philosophy of our 
ansatz supposes that it averages over local cross-section 
fluctuations; and (v) we wish to compare Eq. (jSJ directly 
to experimental data on the charm continuum region such 
as from Beijing |1?tI | . The narrow resonance contribution 



n 


resonances 


continuum 


total 


theory 





1.16 (6) 


-3.03 ±0.37 


-1.86 ±0.37 


input 


( 4) 


1 


1.12 (6) 


1.04 ±0.14 


2.16 ±0.16 


2.19 


( 6) 


2 


1.10 (7) 


0.37 ±0.07 


1.47 ±0.10 


1.49 


( 9) 


3 


1.10 (7) 


0.17 ±0.04 


1.27 ±0.08 


1.26 


(14) 


4 


1.11 (7) 


0.09 ± 0.02 


1.20 ±0.08 


1.16 


(20) 


5 


1.13 (7) 


0.05 ±0.01 


1.18 ±0.08 


1.10 


(31) 





1.17 (5) 


-52.44 ± 1.24 


-51.27 ± 1.24 


input 


( 2) 


1 


1.24 (5) 


3.12 ±0.53 


4.36 ±0.54 


4.51 


( 2) 


2 


1.31 (5) 


1.33 ±0.30 


2.64 ±0.31 


2.79 


( 3) 


3 


1.40 (5) 


0.75 ±0.19 


2.15 ±0.20 


2.27 


( 5) 


4 


1.50 (5) 


0.48 ±0.13 


1.98 ±0.14 


2.06 


( 7) 


5 


1.61 (5) 


0.33 ±0.10 


1.94 ±0.11 


1.99 


(10) 


6 


1.74 (6) 


0.23 ±0.07 


1.98 ±0.09 


1.98 


(14) 


7 


1.89 (6) 


0.17 ±0.05 


2.06 ±0.08 


2.03 


(19) 



TABLE II: Results for the lowest moments, Mn, defined 
in Eq. (0 for n = (t ^ 00) and Eq. ll^J for n > 1. 
The upper (lower) half of the Table corresponds to the 
charm (bottom) quark. Each moment has been multiplied 
by 10"GeV^" (10^"+^GeV^"). The continuum error is from 
AAg'"^ = ±1.47. The last column shows the theoretical pre- 
diction for rhcirhc) = 1.289 GeV, rht{mb) = 4.207 GeV, and 
as{Mz) ~ 0.1211, where the uncertainty is our estimate for 
the truncation error (see text). 



to the various moments is shown in the second column of 
Table im The Ffj are obtained from constrained fits j2| to 
a great number of measurements independently for each 
resonance, and should have very small correlations. We 
will therefore combine their propagated errors in quadra- 
ture. The 3rd column gives the continuum contribution, 
and the 4th column shows the totals to be compared with 
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the theoretical moments in the last column, 



M 



theory 



1 



(6) 



^ 2rhq{faq 

The Cn are known up to 0{a1) and from Refs. 0,0] 
where they were computed for arbitrary renormalization 
scale jjL. It seems appropriate to choose /i = •mq{rhq), 
eliminating all logarithmic terms as there is only one scale 
in the problem. Indeed, the authors of Ref . i who have 
chosen /i = 3 (10) GeV for the charm (bottom) quark and 
then evolved to /i = rhg, report a variation over the first 
5 (7) moments of 122 (312) MeV. (For larger moments the 
as expansion ^3 of the gluon condensate contribution 
breaks down.) Using the same moments |0 but choos- 
ing /i = rhq instead, we observe a variation of less than 
27 (16) MeV. This impressive improvement clearly over- 
compensates for the larger ots- We will choose ^ = mg 
in the following. As for the theoretical uncertainty as- 
sociated with the truncation of the perturbative series, 
we use the method suggested in Ref. jl3| . It exploits the 
fact that once stripped off all group theoretical factors, 
the coefficients appearing in PQCD (to the very least for 
highly inclusive quantities defined in the Euclidean do- 
main) are strictly of order unity (IS]. Since one can easily 
determine the largest group theoretical factor in the next 
uncalculated perturbative order, this offers a reliable and 
transparent way to assess truncation errors. In our case 
this yields the error estimate. 



± NcQ'qCFC. 



1 



2rnq{mq) 



2n 



(7) 



{Nc = Ca — 'iCp — 3) corresponding to ±16ag/7r'^ in 
the Cn- Comparing the corresponding estimate against 
the exactly known coefficients of the first eight moments 
up to order 0, 0| shows that with = m,, 23 of 
24 coefficients are within the estimate, while only one 
coefficient would have been underestimated by a factor 
« 1.437. This seems to be a reasonable state of affairs 
for a 1(7 error estimate and corresponds to ±20 MeV for 
rhrirhf.) from Mi, while variation of the renormalization 
scale [llj assesses this error to only 1 MeV, which is op- 
timistic. We show the estimate {7} in the last column 
of Table m The last two columns of that Table would 
agree within errors even if we had chosen significantly 
smaller variations in A3 and especially A§ (AA3 = ±1.47 
accounts for the error introduced by our ansatz and is 
above and beyond the variations induced by the fit pa- 
rameters). The reason for our more conservative error is 
shown in Table IIIII It shows that Eq. © with A§ = 0.50 
reproduces the n dependence of the moments computed 
from recent data by the BES Collaboration [3 remark- 



^3 



1.97, 



ably well. However, our method favors A§ 
and thus 30 to 40% larger contributions. Table ITTll also 
compares the BES data to the ^'(35') contribution 2] in 
the narrow width approximation. Even assuming that 



n 


BES 


A§ = 0.50 


= 1.97 


BES *(3S') 





5.51 (35) 


5.50 


7.19 


0.215 (39) 0.348 (54) 


1 


3.02 (19) 


3.01 


3.98 


0.151 (27) 0.245 (38) 


2 


1.68 (11) 


1.68 


2.25 


0.106 (19) 0.172 (27) 


3 


0.95 ( 6) 


0.96 


1.29 


0.074 (13) 0.121 (19) 


4 


0.55 ( 4) 


0.55 


0.76 


0.052 ( 9) 0.085 (13) 


5 


0.32 ( 2) 


0.33 


0.45 


0.037 ( 6) 0.060 ( 9) 



TABLE III: The left part shows contributions to the charm 
moments (x lO'^+^GeV^") from 2Moo < < 4.8 GeV, and 
the right part from 2Moo < ^/s < 3.83 GeV. Following 
Ref. [ijl, we computed the columns labeled BES by sub- 
tracting from the threshold data on R{s) the average, R, 
below threshold. (We applied corrections for the leading s- 
dependence.) The errors combine the statistical and uncorre- 
lated systematic ones of R with those in the continuum region 
and with the common systematics (< 3.5%) of the difference. 



the ^'(35') resonance (Af^(35) = 3.7699 GeV) saturates 
the charm cross-section in that region, we observe a di- 
rect experimental 2a discrepancy between Ref. |l5| and 
-'^*(3S) ~ 0.26 ±0.04 keV 1^. Thus there is a discrepency 
between the theory and the BES data, though the the- 
ory does seem to be consistent with the ^(35*) data. This 
constitutes a great puzzle which needs to be resolved in 
the future. We may be able to quote smaller errors after 
this situation has been resolved. Nevertheless, the quark 
masses can still be determined precisely through the sum 
rule approach. 

There is a possible contribution from the gluon con- 
densate J]- It is known up to 0{as) 17], but its actual 
value is not well known. Its inclusion lowers the extracted 
quark masses, increases A§, and sharpens the discrep- 
ancy with the BES data. We can bound its value to 
^ 0.07 GeV^ by demanding n independent results within 
the uncertainties. We use this bound (with a central 
value of zero) to account collectively for non-perturbative 
uncertainties. They induce errors of about 29 MeV into 
Tnc{rfic) {n — 2) and 2.4 MeV into mi,{rhh) {n — 6). 

The parametric uncertainties from as and the quark 
masses themselves are correlated in a complicated way 
(i) across the moments, (ii) across the two quark flavors, 
(iii) between the theoretical moments and the continuum 
contribution, and (iv) with each other. In practice, all 
this is accounted for by performing fits to the moments. 
Heavy quark radiation by light quarks is not res- 
onating and problems associated with singlet contribu- 
tions 0,113] appear only at 0{a^), so these issues should 
not introduce further uncertainties into our analysis. We 
will present our final results after discussing the r life- 
time. 

For our analysis of the r mean lifetime, 



h- 



l-Bs 



290.96 ±0.59 fs. 



(8) 



we evaluate the partial widths into leptons, F^ + r^, and 
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hadrons with vanishing net strangeness, , theoreti- 
cally. The relative fraction of decays with AS = — 1, 
Bs — 0.0286 ± 0.0009 0, is based on experimental data, 
since the value for the strange quark mass, rhs^rrir), is 
not well known, and the PQCD expansion, Cq^^, pro- 
portional to converges poorly and cannot be trusted. 
^QCD multiplies the corresponding to^ ^ terms in 
r"'', posing the same but numerically less important 
problem there. We solved it, by relating to 

the ratio r;i^|Kd| Vlr^'^lKs P) = 0.896 ± 0.034 (in 
which to linear order all universal terms cancel), and 
find C^^liml - ml) = (0.091 ± 0.046). We included 
one-loop electroweak (^ew) [13 and QED (phase space) 
corrections [2^, quark condensate contributions, as well 
as c quark effects in an expansion in m^/4m^ (23 • E.g., 



respectively, as constraints in a fit to all data [3 yields, 

a.(Afz)-0.121ltO:™l?, 

rhrim,) = 1.289tEJ;°l° GeV, (9) 

mb{mb) = 4.207i°:°^5 GeV. 

These results reduce the error in a[Mz) by 25%. 
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10047005. 
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QCD 
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a 85 TT^ 
^*^24 ~ Y' 



4m^±Q!^] 



-pQCD — 1 is the massless QCD correction. Due to ef- 
fects governed by the /3-function, the 3-loop PQCD ex- 
pansion shows slow convergence. The series can be 
reorganized 1241 into a well-behaved expansion with co- 
efficients, di [1^, from the Adler D-function. The new 
expansion is not a power series. Rather, the di multiply 
complicated functions, Ai{as) 24], which we calculate 
numerically up to 4- loop order in the f3 function j25|. 

We computed the world average (jSJ by combining the 
direct value, = 290.6 ± 1.1 fs with r^(6e,6^) ^ 
291.1 ± 0.7 fs derived from the leptonic branching ra- 
tios Be = 0.1784(6) and B^ = 0.1737(6) B\ taking into 
account their 1% correlation. The dominant theoretical 
error induced by the unknown coefficient d^ — 0±77 
is itself strongly aj-dependent, is recalculated in each call 
within a fit, and induces an asymmetric error. 

Other experimental uncertainties arise from 2] = 
1.77699(28) GeV, = 0.97485(46), and Bs. Uncer- 

tainties from higher dimensional terms in the operator 
product expansion, OPE, are taken from Ref. |2y] and 
add up to Ar^(OPE) = ±0.64 fs. We assume that 
an uncertainty of the same size is induced by possi- 
ble OPE breaking effects |2^. The unknown five- loop 
/3-function coefficient, /?4 = ± 579 contributes 
mainly to the evolution of as{mr) to as{Mz) and less 
to the Ai. The subleading errors listed in this para- 
graph amount to ±1.2 fs. We find, asimr) = 0.356lg 
and as{Mz) = 0.122llg;QQ23, in excellent agreement with 
as{Mz) = 0.1200 ± 0.0028 from Z-decays and most 
other recent evaluations of Tt |26l l27j . Including Tr , and 
the n = 2 and n = 6 moments for the c and b quark. 
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